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^ ! Abstract 
OO , 

\ The flow of granular material in a rotating cylinder was simulated by molecular 

dynamics in two dimensions using spherical as well as nonspherical grains. At very 
■ low but constant angular velocity we found that the flow varies irregularly with 

CN , time. The particles move stick-slip like i.e. there are avalanches of different size at 

the surface of the granular material. Observing the traces of the particles we found 
' that there are unstable convection cells. Our results agree with recent experiments 

^ . by Rajchenbach [1] and Rolf [2]. 

o 
o 



1 Introduction 

> ■ 

•i-H 

X 

■ The flow of granular material in a rotating cylinder is of particular interest 

since on the one hand it reveals many interesting to physicists phenomena, on 
the other hand in chemical engineering widely used tumbling ball mills base 
on the motion of iron spheres moving in rotating cylinders. Hence it has been 
intensively experimentally subjected by both engineers [2] and physicists [1,3]. 

The present paper contains results of two dimensional molecular dynamics 
of grains moving in an uniformly rotating cylinder. Recently Rajchenbach [1] 
has shown experimentally that the inclination of sand flowing in a rotating 
cylinder depends nonlinearly on the angular velocity of the cylinder. If the 
angular velocity is smaller than a certain threshold the sand moves stick-slip 
like, i.e. there is no continuous flow but there are avalanches of different size, 
the angle of inclination fluctuates. Using high speed recording of polydisperse 
iron spheres in a tumbling ball mill Rolf [2] observed that the flow forms 
convection cells. They found that the size and the number of the cells depend 
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on the rotation velocity. We will show that our simulations agree with the 
experimental results [1-3]. 

Traditional molecular dynamics simulations of granular material use polydis- 

perse spherical grains, examples can be found in [4-6] and many references 
therein. Applying a recently introduced model for a nonspherical grain [7] 
which deforms elastically under pressure and dissipates energy during colli- 
sions our results agree better with the experiment than the simulations using 
spheres for the case of low angular velocity. Nonspherical grains in molecular 
dynamics simulations of granular material have been introduced already by 
Gallas and Sokolowski [8] in an investigation of the angle of repose of a sand 
heap. In their model a grain consists of two spheres which are connected by a 
stiff bar, their results agree with ours. 



2 Molecular dynamics of granular media 



In the simulations we assume that the grains are ideal spheres. To describe 
the interaction between colliding particles we use a simple Ansatz given by 
Cundall and Strack [5] and Haff and Werner [6]. Two grains i and j at the 
positions and fj will interact only if they touch each other, i.e. if the distance 
between the center points is smaller than the sum of their radii Ri and Rj 

\n-rj\ < Ri + Rj . (1) 



For this case the force between particles i and j moving with the velocities fj 
and fj and rotating with angular velocities and Ctj is given by 

i?,, = j.y.^ + j.|.(0 -1).^ . (2) 

^ \fi-fj\ ^ Vl 0/ \fi-fj\ 



with the normal force 

= y • (i?, + Rj - In - rj\)l + 7iv • mf ■ jrl - fj\ (3) 



and the shear force 

Ff, = min{-75 • mf • \ir,f\ , . (4) 



For the relative velocity 'Cjj' of the particle surfaces one can write 

v^f = (f, - f,) + Ri-ni + Rj ■ (ij . (5) 
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The effective mass of the particles i and j is given by 

pff mi ■ rrij 
mi + rrij 



and the resulting momenta Mj and Mj for the particles i and j are 

M, = • (7) 

and 

M,- = -i^f • i?,. . (8) 

y is the Young modulus, and 75 are the phenomenological normal and 
shear friction coefficients, ^ stands for the Coulomb friction parameter. 

The underlying assumption is that the particles deform each other slightly 
during collisions. To simulate three dimensional behavior we use the Hertzian 
contact force [9] which rises with the power 1.5 of the penetration depth Ri + 
Rj — \Ti — rj\. Eq. (4) takes into account that the particles slide on each other 
for the case that the inequality 

/^•l^yi < \F^,\ (9) 



holds, otherwise they roll. This behavior is due to the Coulomb law. 

— * 

For the integration of Newtons equations of motion due to the forces Fij given 
above we used a sixth order Cear predictor-corrector algorithm [10]. This 
high accuracy is necessary since there is a short range interaction between the 
particles, a small inaccuracy in the calculation of the particle positions may 
cause large faults in the resulting forces. Using algorithms of lower accuracy 
the system would become numerical unstable. The calculation of the rotation 
of the particles is not such sensitive since the angular velocity is very small in 
general, hence we applied a forth order algorithm only. 

In the simulations we use for the parameters Y = lO^kg/s'^, 'Jn = 1-5 -10*^ s~^, 
75 = 3 ■ 10^ s~^, fi = 0.5. The radii of the particles are Causs distributed 
with the mean value Ri = I mm. The chosen parameters are typical for a soft 
granular material. 

The model described above was applied to a simulation of granular material 
consisting of 2400 single grains moving in a clockwise rotating cylinder of 
diameter 160 ■ R^, the cylinder rotates with angular velocity fl = 0.002 sec~^. 
Fig. 1 shows a snapshot of the simulation. 
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Fig. 1. Snapshot of the simulation of 2400 spherical grains moving in a cylinder 
which rotates with constant angular velocity Q = 0.002 sec~^. High velocity is 
coded by red color, particles moving with low velocity are drawn blue. 

The color codes for the velocity of the particles, slow particles are drawn blue, 
fast are drawn red. As visible the velocity of the particles near the surface 
of the flow is nonuniform. The surface of the granular material is close to an 
even plane which is inclined by an angle 0^8°. Although the dynamics of the 
flow seems to be very regular the traces of the grains during the simulation 
are due to convection cells (fig. 2). Since there are more than one cell they 
cannot remain stable over a long time because all of them have the same sense 
of rotation and therefore they disturb each other. This effect was observed 
experimentally too in the investigation of tumbling ball mills [2]. 




Fig. 2. Traces of the motion of granular material at an angular velocity 
Q = 0.002 sec~^. One observes several unstable convection cells. 
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Fig. 3. Snapshot of the simulation of 2400 spherical grains moving in a cylinder 
which rotates with high angular velocity Q = 0.1 sec~^. The bulk of the granular 
material is S-shaped. 



Fig. 4. Traces of the particles for a high angular velocity = 0.1 sec ^). There is 
only one large convection cell. 

Figs. 3 and 4 show a snapshot and the plot of the traces of a simulation with 
angular velocity Q = 0.1 sec~^. The surface of the granular material is not an 
even plane anymore, but it becomes S-shaped. The same shape characteristic 
was observed in experiments by Rajchenbach [1]. In the simulations for high 
angular velocity one finds only one large convection cell (fig. 4) instead of 
several for lower angular velocity, this behavior is due to the experiment [2]. 

After submitting this article we received a preprint by Baumann, Jobs and 
Wolf [11] who applied the interesting simulation method by Visscher and Bol- 
sterli [12] (improved version [13]) to the problem of the rotating cylinder. They 
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Fig. 5. Schematic description of the procedure we used to determine the inchnation 



found convection cells as well. 

Because the surface of the simulated granular flow is not a perfect plane even 
for very slow rotation we have to determine the inclination indirectly. It is 
given by 



where {xcmp , Hemp) are the coordinates of the center of mass point and 
(^cc ) Vcc) is the rotation center of the cylinder (fig. 5). 

As shown in figs. 1 and 3 the profile of the granular material depends on the 
angular velocity. Curve (a) in fig. 6 shows the relation between the angular ve- 
locity VL and the inclination ©. It rises almost linear with the angular velocity. 
Using the Ansatz 



we found for the exponent m ^ 0.8 for the interval < Q < 2 rpm. For larger 
angular velocities the profile is not close to an even plane but S-shaped. 

Although the cylinder rotates with a constant angular velocity the flow of the 
surface particles varies irregularly with time and space, this behavior is shown 
in fig. 7. Each horizontal strip of this figure shows the velocity distribution 
along the surface of the flow. As before red color means high velocity, blue 
means low velocity. The irregular flow was experimentally observed by Ra- 
jchenbach [1] and Briscoe et. al. [3] before. Curves (a) in flg. 9 show the time 
evolution of the inclination and of the particle velocity at the surface, both 
values fluctuate irregularly. 
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Fig. 6. Inclination as a function of the angular velocity O of the cylinder for 
spherical (a) and nonspherical (b) grains. In both cases the inclination rises as 
— 0c ~ r^*^'^, where 9c ~ 8° for spheres and @c ~ 18° for nonspherical particles. 

3 A more complex model 



In the previous section we have shown that the granular flow is irregular, how- 
ever, fig. 9 proves that there are no sharp bounded avalanches. Experimental 
observations, however, show that the flow at the surface of granular material 
in a rotating cylinder is stick-slip like, i.e. there are avalanches of different size 
going down the surface. Our simulations with spherical grains cannot repro- 
duce this behavior. Hence we apply another model where the grains are not 
simple spheres but of more complex shape. 

Each nonspherical particle k consists of four spheres of radius r.^'' located at 
the corners of a square of size L^^^ and one additional sphere in the center of 
the square (fig. 10). The radius r^^ of the middle sphere is chosen to touch 
the others 



[12) 



During each collision between to spheres i and j with radii Ri and Rj which 
might belong to the same grain or to different grains acts a force in normal 
direction 



Y{Ri + Rj - \ fi - fj\) + ^Nrriij^lfi - fj\ 



r,: - Ti 



(13) 



where the effective mass is given by eq. (6). The notation is the same as above. 
In addition each pair of neighboring spheres i and j which both belong to the 
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Fig. 7. (left) Evolution of the velocity of the particles moving at the surface of the 
flow. Each raw is a snapshot. The flow varies irregularly in time and space. 



Fig. 8. (right) Flow at the surface of the granular material in a simulation using 
nonspherical grains. There are successive regions of higher and lower energy due to 
avalanches. The color contrasts are sharper than in fig. 7. 
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Fig. 9. The inclination 0p] and the average velocity v[50 sec~^] of the surface par- 
ticles for spherical (a) and nonspherical grains (b). While v{b) fluctuates drastically 
due to avalanches v{a) is much smoother. 




Fig. 10. The nonspherical particles consist of four spheres at the corners of a square 
and one in the middle of the square which radius is chosen to touch the others. 

same grain feel a force F^-^ due to a damped spring 



= [a{Ci'^ - In - r, I) + 75.f |rl - f 



(14) 



a stands for the spring constant and 7 is the damping coefficient. The term 
C^'^^ is the length of the spring in equilibrium 



a 



(k) 



L^^'> a i, j lie at the same edge 
L('=)/v^ else . 



(15) 



We want to point out here that eq. (4) which introduces static friction into the 
interaction forces is not applied in the simulations using nonspherical grains. 
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Fig. 11. Snapshot of the simulation of 1000 nonspherical grains each consisting 
of 5 spheres moving in a cylinder which rotates with constant angular velocity 
U = 0.002 sec^^. The red drawn grains at the surface show an avalanche. 

We investigated the flow of 1000 nonspherical grains, where the parameters 
are the same as in the previous simulation and a = 10'^ kg/ s^, •ysp = 3- 10"^ s~^, 
lik) i^ik) = 4 and = 3 mm. The cylinder of diameter D ^ 86 ■ L^^^ rotates 
with angular velocity Vt = 0.002 sec~^. 

Fig. 11 shows a snapshot of the simulation. Obviously the inclination is ap- 
proximately twice the inclination observed for spherical grains. The particles 
at the surface move with much higher velocities than the others. This is due 
to an avalanche going down the surface in the moment when the snapshot was 
done. 

The equivalent picture to fig. 7 is shown in fig. 8. The fiow is irregular in time 
and space too but one finds time intervals where the surface particles move 
with high velocities and other intervals where the particle velocities, i.e the 
fiow at the surface, are almost zero. This succession indexes stick-slip motion. 
The time evolution of the average velocity v of the surface particles and of 
the inclination 6 of the surface for nonspherical (curves(b)) as well as for 
spherical grains (curves(a)) is drawn in fig. 9. While in the case of spherical 
grains v varies only slightly for the case of nonspherical particles there are 
drastic fiuctuations of the average velocity due to avalanches of different size. 

For the model of nonspherical particles we investigated the dependence of the 
angle G of the surface on the angular velocity Q too (fig. 6). The measured 
inclinations are approximately twice the inclinations observed in simulations 
with spheres. These values agree better with experimental investigations. For 
the exponent m in — Oc ~ fi'" we found m ^ 0.8. 
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4 Conclusion 



The present paper reports results of two dimensional molecular dynamics sim- 
ulations of granular material moving in a cylinder which rotates with constant 
angular velocity. Using spherical particles and the Ansatz of Cundall and 
Strack [5] for the simulation of static friction we found that the flow fluctu- 
ates irregularly with time. The angle of the material as well as the velocity 
of the grains at the surface fluctuate significantly. For the dependence of the 
inclination © of the granular material on the angular velocity Q, of the cylinder 
we found O — 0c ~ Q^'^. The flow of the particles forms unstable convection 
cells. For low rotation velocity we observed several cells, for high rotation ve- 
locity we observed only one cell. In the latter case the bulk of the granular 
flow reveals a typical S-shape. 

When using nonspherical grains, where each grain consists of four spheres 
located at the corners of a square and a flfth sphere at the center, and where 
neighboring spheres of the grain are connected by a damped spring, we found 
that the flow moves stick-shp hke. That means there are avalanches of different 
size at the surface of the flow. The inclination of the granular material is 
more than two times larger than in the simulations using spheres which agrees 
better with the experiment. 

In experiments was observed that for low angular velocity ft the particles move 
stick-slip like due to avalanches going down the inclined surface, for higher 
velocity one found continuous flow and the characteristic S-shape [1]. The 
inclination depends nonlinearly on the angular velocity. In experiments with 
tumbhng ball mills unstable convection cells have been observed [2] . 

Our simulations reproduced the experimentally measured effects using simple 
assumptions about the behavior of a grain during the collisions. The dynamic 
effects as irregular flow at low rotation velocity, the flnite inclination and the 
existence of convection cells could be reproduced using simple spherical grains. 
To simulate the avalanches at the surface of the flow, which is an effect of the 
static friction between the grains, we applied a more complex model for grains 
which includes nonsphericity. 
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